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An adaption of the convergent close-coupling method 
(CCC) to calculation of differential ionization cross sections 
[I. Bray and D. V. Fursa, Phys. Rev. A 54, 2991 (1996)] 
is analyzed in the context of the Temkin-Poet model. The 
asymptotic scattering wave functions and the unitarity rela- 
tion are given for the the model. Our analysis is used to 
reconcile the recent criticism by G. Bencze and C. Chandler, 
[Phys. Rev. A 59, 3129 (1999)] of the formulation of Bray 
and Fursa. It is concluded the use of "distinguishable" elec- 
trons as proposed by Bray and Fursa is not essential in the 
close-coupling description of ionization. 
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The CCC method was developed initially [|l]|2| to de- 
scribe scattering of electrons from one-electron targets. 
The method made use of Laguerre basis states that are 
are L 2 functions to discretize the target continuum. The 
method provided convergent amplitudes for scattering 
to low-lying discrete states and total-ionization cross- 
sections. The theory is fully symmetrized since it is 
based on a symmetrized expansion of the two-electron 
wave function. Bray and Fursa || used the fact that 
the positive-energy pseudostates were approximations to 
true continuum functions to propose a method for cal- 
culating energy differential cross sections for ionization. 
They did this in a way that preserved "two-particle" uni- 
tarity implicit in the CCC formalism with the L 2 target 
states. This had the consequence of yielding different 
magnitude amplitudes C s (e a ,€b) and C s (eb,e a ) for the 
same physically indistinguishable process of the two ion- 
ized electrons emerging with energies e , e&, depending on 
which electron, a or 6, is represented by a Coulomb wave. 
Bray and Fursa argued that this asymmetry could be in- 
terpreted in terms of "distinguishable" electrons. This 
has been criticized by Bencze and Chandler (BC) who 
argue that the symmetrization property of the ampli- 
tudes is a fundamental tenet of scattering theory and that 
the CCC amplitudes must satisfy this property counter 
to the numerical evidence. They conclude that "The nu- 
merical CCC amplitudes have not, therefore, converged 
to accurate approximations of the exact amplitudes." 

It is the dichotomy between a method of calculating 
ionization processes that seems very impressive in the 
quality of agreement achieved between theory and exper- 
iment, and the criticism of the model on a very funda- 
mental level, that provides the motivation for the present 
work. To make the discussion as transparent as pos- 



sible we choose the Temkin-Poet model (TP) [^|JfJ] of 
electron-hydrogen scattering. This model is the solu- 
tion of the scattering problem assuming spherical aver- 
aging over both electrons and solving for the total an- 
gular momentum zero two-electron wave function. The 
Schrodingcr equation becomes 



K dx 2 



dy 2 



min(x, y) 



+ 2E)V s (x,y)=0, (1) 



with boundary and symmetry conditions 

* s (x,0) = * s (0,y) = 0, 
* s (y,x) = (-l) s * s (x,y). 



(2) 
(3) 



Since Eq. ([l|) is separable in the regions x > y, and x < 
y, the solutions satisfying physical boundary conditions 
may be written as 

oo 

y) = Ui (y)e- ik '* + £ C^u 3 (y)^ x 



+ / de b C (e a ,e b ) l u eb (y)e 
Jq 



-\-ik a x 



x > y. 



(4) 



The solution is extended to the region x <y by using 
the symmetry condition (0). The solution (Q) has been 
written assuming the total energy is above the ioniza- 
tion threshold so the momenta are fixed according to 
E = £i + jk 2 = tj + jk 2 = e a + €b. As usual the momenta 
for the ionized electrons are given by $k*r b \ — e a m- The 



£«(&) 



Uj are the manifold of bound target states of the hydro- 
gen atom target with / = angular momentum. The 
continuum functions are the regular I — Coulomb func- 
tions multipled by a Coulomb phase and normalized to a 
delta function in the energy. It is obvious from the form 
of the solution that the C s coefficients are related to the 
S'-matrix elements since we have written the expansion 
in terms of incoming and outgoing waves. We are not 
concerned here with the problem of solving for these coef- 
ficients; this problem has been addressed in [|||] . Rather 
we wish to establish the form of the asymptotic scattering 
functions and the form of unitarity relation that one ob- 
tains for this model when the ionization channel is open. 
For scattering to discrete states the asymptotic form of 
Equ. (@) is 



3=1 



Cf iUj (y)e+ ik ^. 



(5) 
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The asymptotic form for the ionization requires that 
both particles go to infinity and that each is represented 
by an outgoing wave. Therefore it is necessary to extract 
these forms from the integral term in (^). Normalizing 
the u~ b to a delta function in energy asymptotically 



u~ (y) y-oo J -^-e~ lS " sin(fc h y + — In 2k b y + S c ). (6) 

The integral can be readily evaluated as x, y — > oo using 
stationary phase methods (a very readable exposition is 
given in ref. JtJ ) . It is determined by the integrand in the 
vicinity of the stationary points (see for example 
of the terms e lk ^ x±lk >>v _ One finds there is no stationary 
point for the exponential with ~ik b y, and hence no con- 
tribution to the asymptotic form. The only contribution 
comes, as expected, from the outgoing wave giving 



*L(x,y) ~ C s {e ai e b ) l e^l\2E)^ cosar" 1 ^ 



xe y k t, by> , x > y. 



(7) 



Here r = (x +y ) 2 and the energies e a .b are fixed by the 
stationary point according to 



tan a = — = — . < a < — . 

x k„ 4 



The discrete scattering sector corresponds to a = 0. 

We now derive the unitarity relation using the conser- 
vation of flux, incoming and outgoing, through a closed 
surface outside which the asymptotic forms for the wave 
function are appropriate. The current for the TP model 
is defined as j(x, y) = S^ 5 * V'S 5 . The surface we choose 
to apply the flux conservation law to is shown in Fig.l 
and is formed by x = 0, y = 0, x = a and y = a. a is 
assumed to be large enough to ensure use of the asymp- 
totic form for ionization is justified. It is also useful to 
define two strip regions running along the x— and y— 
axes bounded by x — a n , y — a n where \ < r/ < 1. 
It is assumed that the flux for the discrete target state 
wave functions is contained within these strips for a fi- 
nite number of the states which we may choose to be 
as large as we wish. Conservation of flux requires that 
J n j.dn = 0. On account of the boundary conditions (||) 
and symmetry condition (||) this reduces to 



jx(x = a,y)dy = 0. 



(9) 



The incoming incident and the discrete outgoing currents 
are confined to the strips and as such we see that 



/(mc) = -hi \ui(y)\dy ->• -hi, as a -> oo, (10) 
Jo 

oo 

f(disc) ^+Y1 k i\ C 3i I" as a °°- ( U ) 

Calculation of the ionization flux requires more care be- 
cause the asymptotic form of the Coulomb wave m) is 



valid only when k b y >> 1. Thus in order to calculate 
the ionization flux we have to ensure that this condition 
is met throughout the interval of integration [a v , a] . (We 
note that there is a minimum value of kb fixed by the 
stationary phase condition (||).) These constraints may 
be satisfied by choosing 



~ <r?<l, a»(2E)-i. 



(12) 



The total ionization flux using the current from Eq. (0), 
is in the limit of large a 



(225) * 



C 5 (e a ,e 6 ),| 2 — 
ai [a 



y 



-dy. (13) 



The remaining terms in the current involve overlaps be- 
tween ionization and discrete wave functions and are 
readily shown to contribute vanishing flux in the limit 
of large a. Collecting all the flux contributions the re- 
quired unitarity relation in the limit a —* oo is 



1 = 



3=1 



|C s (e a ,e h ) l p 



kb ki 



-de h 



(14) 



(8) Thus the ionization contribution to unitarity requires 
C s (e a ,e b )t only for < e b < f . 

Let us now turn to the close-coupling (CC) approach 
and its description of ionization. One begins with an 
explicitly symmetrized wave function: 



^ cc \x,y) = (l + (-l) s P(x,y))(J2fU x ) u M 

3=1 

de b f s (e a ,e b )i(x)u ei (y)) (15) 

where P is the permutation operator for x, y and 

ffi x ~oo S lje ~^ x + cf/ cc V^, (16) 

(17) 



The asymptotic form of the CC wave function can be 
found using the methods already discussed. It is readily 
seen that the discrete scattering boundary condition (|5j) 
is unaltered but the ionization asymptotic form becomes 



^ S ion C \ x >y) - ^/\2E)i cos ar^/V 



'2Er 



+ {-l) S e^ la2k « x C s{cc \eb,e a )). (18) 

This form has the unsatisfactory feature that it cannot 
be expressed in the form of an outgoing wave multipled 
by an amplitude that depends only on e a , due to the 
logarithmic phases. One is forced to conclude that the 
symmetrization, on its own, does not lead to a generally 
acceptable form for ionization. If the long-range loga- 
rithmic phase were not present, the asymptotic form is 
perfectly acceptable since one may write 
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C S(CC) { 
sym \ 



e b ) = C s ^ cc He a ,e b ) + (-1) S C S < CC ) (e b , e a ) (19) 



and C\ 



S(CC) 

sym 



(e a ,£b) = (-l) S C S y m (e b ,e a ) as demanded 



by conventional quantum theory for identical particles as 
emphasised by Bencze and Chandler pL Further, note 
also that there is no unique form for the f s ^ cc '(e a , e b ) in 
Eq. ©, and hence for the C s< - cc > (e , e&) of Eq. @. 
Thus the iV — ► oo limit in BC cannot be taken as unam- 
biguously as they assume and the CCC equations do not 
have to converge necessarily to a symmetrized ionization 
amplitude. 

Let us now examine how to apply these considerations 
to the CCC method for ionization. In this method the ex- 
pansion over the complete set of target states is replaced 
by a set of L 2 target states generated by diagonalizing 
the Hamiltonian in an TV-function subspace of the La- 
guerre basis. The limit N — + oo gives the CC expansion 
(|l5|) . The asymptotic form of the wave function is 



S(N) 
disc 



(x,y) x - 00 ul N \y)e- ik i x + 

N 
3=1 



(20) 



The Uj are the L 2 approximate target states and their 



Because the vS are L 2 states the 



(N) 

energies are . _ 

scattering flux for all the states is confined to the strip 
regions in Fig 1. Thus as it stands there is no flux in the 
ionization sector < a < ^ . In order to sensibly extract 
an ionization amplitude from Eq. (^0|) one must make 
an approximation that yields a true continuum function 
for the positive-energy target states. The following is an 
obvious possibility: insert the completeness relation for 
the true target states in the variable y and let y — > oo. 
Then one may write 



^\x,y) M ^oo J de b (1 + (-l) s P(x,y)) 

j:0<ef> <E 



where use was made of the fact that (u eb \Uj ) is highly 



(21) 



peaked about ej^ and energy conservation to replace 
e + lk \ >x |-,y e +tk a x_ is straightforward to prove this 
leads to the expression of Bray and Fursa ||] : 

C S {\k^\e^) i ^C^\u^\uf ) )- (22) 

These amplitudes may be interpolated to provide the 
(unsymmetric) C s ^ N ^ cc \e a , e b ) which are then used to 
construct the symmetrized amplitude (|l9|). It is this 
amp litude that must be used in the unitarity relation 
( |i"4| ) . Bray and Fursa |3| one the other hand argued that 
the unsymmetric amplitudes must be summed incoher- 
ently to satisfy unitarity, a suggestion that is equivalent 



to treating the two electrons as distinguishable and the 
integration is performed over the interval < e b < E. 
Bray further suggested jl0[ that in the limit N — > oo the 
CCC(N) amplitudes for? < e b < E would converge to 
zero. This has also been criticized by BC and rebutted 
by Bray fn]] . We now give a specific example where the 
distinguishable-electron hypothesis can be compared di- 
rectly with coherent summation of identical-particle am- 
plitudes. 

Consider the letter of Bray @ that reports CCC(N) 
calculations for the TP singlet scattering model at 
E=2Ryd. It is stated there that " there is a clear lack 
of convergence in going from 10 to 50 states". However, 
curiously, on inspection of his Fig. 2 |l(| there appears 
to be some convergence in the region of y. We can use 
the fact that the CCC(N) cross sections are convergent 
there to extract a sequence of interpolated CCC(N) val- 
ues for the -j cross section. In Table 1 we show the re- 
sults obtained by a spline interpolation. The convergence 
is apparent. The reason this energy is critical is that it 
is the one energy where the symmetrized CC amplitude 
(13) properly satisfies the boundary conditions of Eq. 
(W. Thus the symmetrized amplitude leads directly to 
a cross section proportional to 4|C 0C ' c ' c '( Ar )(-|, f )| 2 - In 
the table are also shown the results of a very recent cal- 
culation by Baertschy Rescigno and McCurdy E^] using 
an exterior complex scaling method 13 that circumvents 
the problem of matching to Coulomb ionization bound- 
ary conditions. The result through use of the correctly 
symmetrized CC amplitude is in good agreement with 
their value. The symmetrized CC amplitude of course 
is symmetric about -j and no step occurs at -j- As we 
move away from this energy on the assumption that the 
neglect of logarithmic phase factors is a secondary effect, 
the cross section should move smoothly up from its min- 
imum value. This is the behavior observed by Baertschy 
ct al jl2|. 

Interestingly the the total ionization cross section ob- 
tained by Baertschy et al is in excellent agreement with 
the CCC predictions |Q. Using the symmetrized am- 
plitude and integrating up to seems to lead to a very 
similar answer as using incoherent summation and inte- 
grating from < e b < E. Thus all the features of the TP 
model are explained by CCC theory provided one uses 
a fully symmetrized ionization amplitude. In case it is 
argued that the example chosen is fortuitous to our case, 
we also analyzed the TP cross-section data presented by 
Bartschat and Bray |1|| at E=3Ryd. The results are 
shown in Table 1. Again there is convergence to the 
value of Baertschy et al. 

Our conclusions are as follows. For CC methods in 
general, the symmetrized expansion gives rise to ion- 
ization scattering wave function asymptotic forms that 
cannot be simply expressed as outgoing wave in hyper- 
spherical coordinates multiplied by an amplitude unless 
there are special circumstances. For TP triplet scatter- 
ing one such special case arises in that CCC calculations 
yield an unsymmetric amplitude that converges to zero 
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for £{, > -j. In this case a symmetrized amplitude can 
be made trivially. Thus if the TP CCC triplet scattering 
calculations converge they must converge to the correct 
model values and this is what is observed For TP 
singlet scattering at the equal energy -sharing kinemat- 
ics, efc = e a = -j, the symmetrized CC amplitude satisfies 
the correct ionization asymptotic form including loga- 
rithmic phases. Using the converged unsymmetric sin- 
glet CCC amplitudes to construct the symmetrized CC 
cross section we get very good agreement with an inde- 
pendent method of solution for the Temkin-Poet model. 
These kinematics are ones where coherence of amplitudes 
is maximal. The correct symmetrized singlet cross sec- 
tion is therefore a factor of four larger than the value 
from the raw CCC calculations. For asymmetric energy 
kinematics the symmetrized amplitude will be very sim- 
ilar to the published CCC calculations below -j because 
the unsymmetric amplitudes above -j are small. 

We make the following observation regarding Bray's 
step hypothesis. It is classical Fourier theory behaviour 
that Fourier series which represent functions with steps, 
converge, at the point of discontinuity, to the midpoint 
of the step. If the singlet CCC amplitude has such a 
step, one might expect the amplitude to converge to the 
midpoint at the step and the cross-section therefore to 
one quarter of its step height. This is consistent with our 
findings. 

The conclusions regarding coherent summation of am- 
plitudes may carry over to the full e-H and e-He prob- 
lems. For example, it has been found that the CCC 
results, after incoherent combination, at equal-energy- 
sharing are typically a factor of two too low, i.e. a fac- 
tor of four on the raw CCC results Jl|. For asymmetric 
kinematics it is likely that the CCC converges to the true 
amplitudes contrary to the analysis of BC. 

The author is indebted to Igor Bray, for stimulating 
discussions and points of clarification regarding his work 
and for providing access to his data base of CCC calcu- 
lations. The support of the Australian Research Council 
and Murdoch University is also acknowledged. 
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TABLE I. Singlet energy DCS (in units of 7ra„ per Ryd) 
of the TP e-H scattering problem are shown at a total 
three-body energy of 2Ryd (To) and 3Ryd [[jj for CCC(N) 
calculations with the kinematics = e a = ^. They are mul- 
tipled by a factor of 4 for reasons given in the text. N pos is 
the number of positive energy open-channel pseudostates in 
each calculation. The CCC(N) numbers are compared with 
the calculation of Baertschy, Rescigno and McCurdy |12[. 
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(20,8) 
0.0172 


(30,13) (40,18) 
0.0145 0.0137 


(50,24) 
0.0135 


0.0140 
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(42,19) 
0.0072 


(44,20) (46,22) 
0.0081 0.0069 


(48,23) 
0.0069 


0.0068 




FIG. 1. The figure shows the closed rectangular sur- 
face and incoming, discrete and ionization currents that pass 
through the surface. See the text for details. 



4 



